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Abstract 

Modem parton distribution function fits require an accurate treatment of heavy quark mass effects. In this contri¬ 
bution, we discuss how the FONLL general-mass variable flavour number scheme can be extended to allow for the 
possibility of a fitted charm PDF. We present a first estimate of the impact of a fitted charm contribution in the FONLL 
structure functions, and we discuss the development towards including a fitted charm PDF in the NNPDF framework. 
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1. Introduction 

One of the main sources of systematic errors in 
hadron colliders such as LHC is the uncertainty that 
affects the parton distribution functions (PDFs), which 
limit the accuracy of our knowledge of proton structure. 
Given the increasing precision of experimental data, it 
is necessary to have a reliable estimate of the PDF un¬ 
certainties. 

Currently, three major collaborations provide global 
PDF sets. Their most recent releases are NNPDF3.0 IfTI . 
MMHT14 121 and CT14 ID. A detailed comparison was 
performed in [51 and shows a reduced discrepancy 
between these sets with respect to previous releases. 

A potential source for the differences in current PDF 
fits is the treatment of heavy quarks. Modern PDF sets 
include heavy quark effects by resorting to general-mass 
variable flavour numbers (GM-VFN) schemes. In cur¬ 
rent PDF fits the GM-VFN schemes are used under the 
assumption that the heavy quark is generated perturba- 
tively above the threshold. Whereas this can be a rea¬ 
sonable assumption for top and bottom, one might ques¬ 
tion whether it is possible to rule out the existence of 
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non-vanishing charm PDF below threshold, sometimes 
called intrinsic charm (IC). Moreover, even if the charm 
PDF is entirely perturbative, and therefore vanishes be¬ 
low the threshold, the computation of the massive cor¬ 
rections in any GM-VFN scheme shows a substantial 
dependence on the charm mass due to the large value 
of asinic). The inclusion of a fitted charm PDF provides 
a solution for both these problems. In this contribu¬ 
tion, we will present some preliminary results showing 
the structure functions for neutral current deep inelas¬ 
tic scattering (DIS) obtained using an extension of the 
FONLL GM-VFN scheme ||6l|71 which includes a fitted 
charm contribution. 

2. NNPDF3.0 

NNPDF3.0 is the latest PDF set released by the 
NNPDF collaboration. The PDF set features new data, 
improved methodology, and a new analysis code. 

Besides the data included in previous releases, sev¬ 
eral new data are added to the dataset. These include 
HERA-II deep-inelastic inclusive cross-sections from 
ZEUS and HI, and the HERA combined charmed pro¬ 
duction cross-section data. NNPDE3.0 moreover in¬ 
cludes a large amount of new data from LHC such as jet 
production from CMS and LHCb, top quark pair pro¬ 
duction total cross section from ATLAS and CMS, and 
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W + c data from CMS, which provide important infor¬ 
mation about strange PDFs. 

NNPDF3.0 features many improvements in the fitting 
methodology. In order to include and process the new 
data, the fitting code has been completely rewritten in 
C-n- optimising the computational intensive hadronic 
calculations used in the fitting methodology. Addi¬ 
tionally, the minimisation of the genetic algorithm uses 
a new mutation approach, both improving the fitting 
speed and the fitted results. NNPDF3.0 also presents 
an extended kinematic range for positivity observables, 
thus obtaining better constraints. Finally, NNPDF3.0 
introduces an improved form of cross-validation, which 
prevents over-learning in the neural networks while also 
reducing the chance of under-learning due to a prema¬ 
ture stopping of the fit. 

The benchmarking of fitting methodology becomes 
particularly important due to the increase in the exper¬ 
imental data included in NNPDF3.0 and the increased 
precision the resulting PDFs achieve. As data becomes 
more abundant, covering a wider kinematic range with 
higher precision, it becomes essential to eliminate any 
uncertainties due to the methodology, such that the 
PDFs reflect theoretical and experimental uncertainties 
exclusively. 

NNPDF3.0 uses closure testing for methodological 
validation. The idea behind closure tests is rather sim¬ 
ple: take an assumed form of the PDFs, a given theo¬ 
retical model (pQCD), and then generate a set of global 
pseudo-data with realistic statistical properties. If the 
fitting methodology is performed using the pseudo-data, 
this should reproduce the assumed form of the underly¬ 
ing PDF, within the correct uncertainties. Closure test¬ 
ing allows one to assess directly the methodology; once 
validated, the methodology can be used in real fits. 

All the NNPDF3.0 fits are available on LHAPDF, 
with results based on LO, NLO and NNLO QCD the¬ 
ory for various values of a^. 

3. Fitted charm 

There are two main motivations for a fitted charm 
PDF. The first is related to the current treatment of the 
heavy quark PDFs in global fits. In existing published 
parton fits the heavy quark PDF is generated dynami¬ 
cally above the threshold. Whereas the bottom and top 
quarks are sufficiently heavy and thus can be treated 
entirely perturbatively, the charm sits at the border be¬ 
tween perturbative and non-peiturbative behaviour and 
in principle one should account for a non-perturbative 
contribution to the charm PDF close to the threshold. 
The second motivation is that even if most of the charm 


PDF is of perturbative origin, it is sensitively dependent 
on the choice of the scale at which PDFs are param- 
eterised. The inclusion of a fitted charm PDF would 
therefore ensure unbiased results since the ignorance 
about the vanishing scale would be absorbed by the 
PDF. 

Over the years, several non-perturbative models of in¬ 
trinsic charm have been proposed iHIol (see e.g. HD 
for a recent review). Currently, the general belief is 
that the intrinsic contribution is of the order of 1% or 
less and thus can be ignored in the majority of collider 
applications.. Charm present at the initial scale may 
have however some very distinct signatures in LHC pro¬ 
cesses, such as direct photon in association with charm 
jets ini or open charm production ISlITl. 

In a recent study ca the CTEQ collaboration per¬ 
formed a global fit that incorporates different IC models 
into the charm PDF. In future releases, NNPDF will in¬ 
clude a fitted charm distribution that will however not 
depend on any IC model, thus ensuring an unbiased re¬ 
sult. 

In the NNPDF framework, heavy quark effects are 
taken into account using the FONLL GM-VFN scheme. 
The scheme provides a consistent matching between the 
fixed order (or massive) calculation, in which the heavy 
quark is treated as massive, and the resummed (or mass¬ 
less) calculation, where the heavy quark is treated as 
massless and the collinear divergences are factorised in 
the PDFs. In I?) three variants of the FONLL scheme 
were proposed. Each variant is defined by the order of 
the resummed calculation and the number of terms in its 
expansion replaced by their massive counterparts. In the 
EONLL-A scheme one considers the resummed calcu¬ 
lation at NLO with one term replaced; in the EONLL- 
B scheme the resummed calculation at NLO and two 
terms replaced; in the EONLL-C the resummed calcu¬ 
lation at NNLO and two terms replaced. 

In order to include a fitted charm, the formulation of 
EONLL as formulated in Q must be modified by extra 
contributions, which are formally subleading in the case 
the charm PDF vanishes below the threshold. While the 
old implementation was equivalent to S-ACOT Ha [13, 
the new one is equivalent to ACOT IfTSl . In particular, 
it is necessary to include LO and NLO diagrams with 
a massive quark in the initial state to take into account 
the possibility of an initial state charm. Such diagrams 
incorporate the mass effects with incoming heavy quark 
lines. These corrections were first computed in lfT9l and 
then extended to general masses and couplings in ||20l, 
which also corrected some errors present in nil. We 
have confirmed that the results of ll20l properly correct 
the original calculation of im. 
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The extension of the FONLL scheme to include fitted 
charm will be discussed in a forthcoming paper ETIl . 
Here we present some preliminary results by consider¬ 
ing the charm structure function and showing the ef¬ 
fect due to the additional terms that must be added to 
the FONLL expression in the presence of a fitted charm 
PDF. 

In order to quantify the possible effects due to a 
non-vanishing IC PDF we consider a toy PDF with an 
intrinsic charm component as predicted by the BHPS 
model ||8l . 

c(x) = Ax^[6x(l + x) ln(x) -l- (1 - x)(l -l- lOx 4- x^)], (1) 

where A is fixed by the magnitude of the intrinsic contri¬ 
bution. We impose the constraint such that the momen¬ 
tum fraction of the intrinsic charm component is 0.5% 
of the total momentum. The toy PDF was produced 
starting with a NNPDF3.0 NLO PDF with as - 0.118 
using APFEL ll22l . which we modified by adding the 
new matching conditions at the threshold /ic - = 

1.275 GeV. The inclusion of the additional terms to 
standard FONLL is performed with a stand-alone code. 

In Fig. we show the structure function in sev¬ 
eral schemes as a function of the energy scale for 
different values of x. On the left we show the results ob¬ 
tained with NNPDF3.0 NLO, whereas on the right we 
show the results obtained using the toy PDF, which in¬ 
cludes an intrinsic charm component. For each value of 
X, we compute F^ in the massive scheme, in the mass¬ 
less scheme, in the standard FONLL-A scheme and in 
the extended FONLL-A scheme. 

We can observe that the effect of the additional 
terms that were added to the FONLL implementation of 
Ref. Q in order to make it compatible with fitted charm 
correspond to a tiny correction to the original FONLL 
result Q at all the values of x considered in the case 
of a vanishing IC PDF. In presence of an intrinsic com¬ 
ponent of the charm PDF, the impact of the additional 
contribution is not negligible anymore for 2 < 10 GeV; 
since the BHPS model Eq. Q that we considered has a 
peak at X ~ 0.2 this is particularly evident at large values 
of X. Eor smaller values of x, where the IC contribution 
in the BHPS model is smaller, the effect of the addi¬ 
tional terms is minor. In this case, the non-vanishing 
scenario becomes closer to the vanishing scenario, as 
it can be seen by comparing the plots in the two cases 
at X = 0.1, X = 0.05. However, the pattern described 
above is due to the fact that we assumed a particular 
model for the intrinsic charm PDF Generally, we ex¬ 
pect to find a different result when the charm contribu¬ 
tion will be fitted to the data. 


4. Outlook 

In order to quantify precisely procedural uncertain¬ 
ties in the determination of the PDFs, NNPDF3.0 intro¬ 
duced important methodological improvements. Hav¬ 
ing these uncertainties under control, it becomes of in¬ 
creasing importance to assess other possible sources of 
uncertainties. 

The uncertainty about the hypothetical scale in which 
the charm PDF vanishes makes the assumption that the 
charm PDFs is generated perturbatively at the threshold 
a potential source of bias, which would be solved by 
introducing a fitted charm PDF. 

A global analysis in the NNPDF framework includ¬ 
ing a fitted charm PDF is on the way. Thanks to the 
flexibility and robustness of the NNPDF methodology, 
it will be possible to include a model-independent fitted 
charm PDF in future fits. 
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Figure 1: The charm structure function calculated in different schemes as a function of the energy scale Q, for three different values of x. The 
plots on the left are obtained using NNPDF3.0 NLO where the charm PDF is generated dynamically at the threshold. The plots on the right are 
obtained assuming an IC component given by Eq. 
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